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Abstract 



< 

In this note we analyse the propagation of a small density perturbation in a one-dimensional com- 

i-G pressible fluid by means of fractional calculus modelling, replacing thus the ordinary time derivative 

j^ with the Caputo fractional derivative in the constitutive equations. By doing so, we embrace a vast 

C^ phenomenology, including subdiffusive, superdiffusive and also memoryless processes like classical 

I " I diffusions. From a mathematical point of view, we study systems of coupled fractional equations, 

leading to fractional diffusion equations or to equations with sequential fractional derivatives. In this 

T-H framework we also propose a method to solve partial differential equations with sequential fractional 

^ derivatives by analysing the corresponding coupled system of equations. 

I -J Keywords: fractional calculus modelling, systems of fractional equations, fractional diffusion-wave 

(<— N equation, fractional sequential derivatives. 

^ 1 Introduction 

en 

,__! Fractional calculus modelling has recently had different developments in the applied science II13II . 

;;• The Caputo derivative of real order introduces a memory formalism as it is an integro-differential 

. ^H operator defined by the convolution of the ordinary derivative with a power law kernel. In this note we 

^ propose an analysis of the propagation of a small density perturbation through a compressible fluid 

J-j using techniques related to fractional calculus. In practice we replace the ordinary time derivative with 

the Caputo fractional derivative in the constitutive equations aiming at furnishing the model with a 

memory This approach leads to a coupled system of fractional differential equations that has to be 

properly treated. 

From a physical point of view, this means introducing a memory formalism in a simple model of 
sound propagation through a one-dimensional compressible fluid. By means of this model we can take 
into account a time delay in the wave propagation. Similar fractional generalisations in modelling 
sound propagation were already considered in some papers. For instance in some recent works, Fellah 
et al. l3j ^ studied transient-wave propagation in porous materials using fractional modelling to take 
into account the frequency variability of some dynamic cofficients of the medium like tortuosity and 
compressibility. Their approach is also validated by experimental observations O . 



Furthermore, we have also to notice that Tarasov [14] gave a space-fractional formulation of the 
hydrodynamic equations to describe fluid flow in fractal media. He also studied sound wave propagation 
within this conceptual framework. In our formulation we use a time-fractional generalisation of the 
characteristic master equations; therefore our model is in some way similar to the one of Fellah et 
al. m concerning propagation in porous media, but more general, as we are going to discuss. As a 
matter of fact we want to show that by using our approach, with some specific hypotheses on the range 
of variability of the fractional order of differentiation, a fractional diffusion-wave equation will result 
by decoupling the constitutive equations. This result appears interesting because we obtain a family 
of equations, depending on the real order of differentiation, which includes the heat equation as a 
particular case. We find thus that, by introducing fractional derivatives in the constitutive equations 
of a one-dimensional compressible fluid, we embrace a vast phenomenology, including subdiffusive 
and superdiffusive processes. Moreover, with our mathematical hypotesis, processes like classical 
diffusions, can be considered as special cases of a more general and complex phenomenology that can 
be treated with fractional calculus. This means that in some physical problems, a diffusive equation 
corresponds to a microscopic fractional process of migrating fluid particles with memory. From a 
mathematical point of view, we study different cases: when the semigroup property is valid for the 
Caputo fractional derivative, we derive a fractional diffusion equation by decoupling the related system 
of equations; otherwise we obtain partial differential equations with fractional sequential derivatives. 
The organisation of this paper is the following: in Section|2]we recall some basic concepts and properties 
related to fractional calculus; in Section [s] our model of propagation of a small perturbation with time 
fractional derivative is discussed; in Section |4] we provide a mathematical discussion on the range 
of variation of the characteristic indices of differentiation; Section [s] concerns the solution of initial 
value problems involving fractional sequential derivatives. This is particularly important as they arise 
naturally for some values of the characteristic indices of fractionality; finally, in Section [6] we give a 
complete discussion of the results in relation to the previous research on generalised sound propagation 
models. 

2 Basic concepts and theorems 

In this section we recall some definitions and basic results related to fractional calculus. 
Definition 2.1. Let a e R"*". The Riemann-Liouville fractional integral is defined as 

"^r/(o=-^ (t-Tr-i/wdT, t>o. (2.1) 



Note that this operator is such that J°/(f ) = /(f )• Moreover, it satisfies the semigroup property, i.e. 

rjpf[t) = r+i^f[t), a > 0, /3 > 0. 

Different definitions of fractional derivatives have been introduced and studied in the literature 
(see for example [11] ot [8]). For the sake of our analysis and because it proved to be of practical 
application, we make use of the Caputo fractional derivative. 

Definition 2.2. Let m — l<a<m, with m e N, the Caputo fractional derivative is defined by 



\fr-'"Kt), 



TTr T „ V- ~ !■ J / I t JU t , U. -f- III, 

D"/(t)=J"'-"D"'/(t)={ r('T.-a)Jo^ J y J ' ^ ' ^2.2) 



where f^^'it) = ,„ is the ordinary derivative of integer order m. 



It is clear by Definition 2.2 that the fractional derivative is a pseudo-differential operator given by 



the convolution of the ordinary derivative with a power law kernel. Therefore the reason why fractional 



derivatives introduce a memory formalism becomes evidem. It is simple to prove the following properties 
of fractional derivatives and integrals (see e.g. Illlin : 

D«J«/(t)=/(0, a>0, (2.3) 

m-l k 

j;^D«/(t)=/(0-2/W(0)-, a>0,t>0, (2.4) 



k=0 



r(fi + 1) „ 

DftP = —^ ' tP-" a>0, /3eR, t>0. (2.5) 

Finally, a key point for the following discussion is the limit of validity of the semigroup property 
for the Caputo fractional derivative. In its recent book, Diethelm (see [2j, pag.56) gave a sufficient 
condition for the validity of the semigroup property for the Caputo fractional derivative. Let us recall it. 

Theorem 2.1 (Law of exponents). Let f(^t) & C'^[0,x], for some x > and some k &N. Moreover, let 
a,l3>0be such that there exists some I e N with I < k and a, a + ^ e [Z — 1, Z]. Then 

D^D'^fit) = Dt^Pfit). (2.6) 

The proof of this theorem is very simple. It proceeds by considering the definition of the Caputo 
derivative in relation to the Riemann-Liouville fractional integral. This theorem highlights a constraint 
on the applicability of the semigroup both with respect to the request of smoothness of the function and 
with respect to the ranges of the real orders of differentiation a and p. This means, for example, that, if 
a G (0, 1], then the law of exponents is applicable if ^ e [0, 1 — a) and /(t) e C^, with fc > 1. Here 
we have also to notice that in most cases the law of exponents is not applicable for fractional Caputo 
derivatives, but anyhow there are different techniques to handle sequential fractional derivatives (see 
for example [lll|). It is clear that in cases where the semigroup is still valid for the Caputo derivative, 
the fractional equations are more manageable and, what is most important, the construction of a Cauchy 
problem demands conditions with meaningful physical sense. In the following discussion therefore, 
we develop our analysis in the cases in which the semigroup property is valid, neglecting the analysis 
of the fractional sequential case for which the requested initial conditions have not a clear physical 
interpretation. 

3 Propagation of a small perturbation by means of fractional cal- 
culus modelling 

Let us consider a one-dimensional compressible fluid, initially at rest, with an homogeneous constant 
density Pq. We study the effect of a small perturbation in the density pXz, t). We thus deal with a 
perturbative theory with 

p(z,t) = po-l-p'(z,t), wliz,t) = wXz,t), zeR, t>0, (3.1) 

where w[z, t) is the velocity field and where p'(,z, and wXz, t) are first order terms. This is a simple 
way to obtain the sound wave propagation equation in a compressible fluid. We will develop this 
classical model replacing the ordinary derivative with the fractional Caputo derivative. 
The equations of continuity and velocity of the fluid, linearised at the first order, are 



d^'pXz, t) + Po5,w'(z, t) = 0, a > 0, /3 > 0, t > 0, 
PodfwXz,t) + cfd,pXz,t) = 0, 



(3.2) 



where cf = dp' /dp' is the sound velocity in the fluid. This is a coupled linear system of fractional 
equations. Note that we have used the symbol d" for the Caputo firactional partial derivative and in 



general the real order a and /3 can be different. In the following we discuss the range of variability of 
the fractional orders in R. 

Now, applying the fractional derivative of order ji to the first equation, and differentiating the 
second equation with respect to z, we obtain 

id^d^pXz, t) + Po3,3f w'(z, t) = 0, 
Wdfw'iz, t)+ i3,,p'(z, t) = 0, 

V Po 

and finally 

3f3>'(z,t)-c2a,,p'(z,0 = 0, a>0,p>0,t>0. (3.4) 

In general this is a partial differential equation with a time fractional sequential derivative (see 111111 1. 
We are interested in the cases in which the law of exponents (Theorem |2.1[ l is satisfied and the equation 
(13. 4D therefore becomes 



S^+'^p'Cz, t) - cfd,,pXz, t) = 0, t > 0, z e R, (3.5) 

i.e. a fractional diffusion-wave equation (when a + ^ e (0, 2]), extensively studied in literature (see 
for example f9l for the fundamental solution) . In fact we wish to arrive at a meaningful equation that 
can be treated considering initial conditions with physical meaning, differently from the case usually 
treated with fractional sequential derivatives. In the next section we thus discuss the cases in which this 
coupled system leads to a single fractional diffusion-wave equation. Note also that with a symmetric 
manipulation we find the following equation for the velocity field wXz, t): 

dyfw'liz,t)-cfd,^wXz,t) = 0, t>0, zeR. (3.6) 

4 Coupled systems of fractional equations related to diffusion and 
subdiffusion 

In the previous section, we showed that our coupled system of fractional equations on velocity and 
density fields of a one-dimensional compressible fluid, is related to the fractional diffusion-wave 
equation, under the hypothesis that the semigroup property holds. Here we discuss in detail the cases 
in which we can apply the law of exponents. This section is devoted to mathematical analysis; we will 
return to physical meaning in the last section. 

We first consider equation ( |3.4D with the constraints 

ae(0,l), 13 & [0,1 -a], (4.1) 

i.e. the pair of parameters (a, P) varies in region ^4 shown in Fig.[l] Note that the constraints ( |4.1| ) are 
equivalent to 

/3e(0,l), ae(0,l-^]. (4.2) 

When ( |4.1| ) or ( |4.2| ) hold, for the Caputo fractional derivative a semigroup rule is valid (see Theorem 
|2.1||, then we can rewrite (13.4|) as 



?«+/3„'r^ ^^_.2;, „/r„ ^^_n ""tL^ ^r„u„ ^^_^2' 



d^' pXz,t)-cfd,,pXz,t) = ^ 5/p'(z,t)-c25,,p'(z, = 0. (4.3) 



Equation ( |4.3| ) is clearly a time-fractional diffusion equation (see [llj, page 296). Conditions ( |4.1| ) 
implies that j e (0, 1], therefore the regime is subdiffusive or at most diffusive. The diffusive regime 
(corresponding to the classical heat equation) is reached for each pair (a, P) for which p = 1 — a (see 
Fig.[l]). 
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We are now going to study the possibility of arriving at a fractional diffusion equation displaying a 
superdiffusive behaviour. In this case, in order to maintain at least partially the semigroup property, we 
consider the following constraints for the fractional indices a and /3 : 

ae(l,2), ^e(0,2-a]. (4.4) 

In Fig. [T] the region D of the a/3 -plane in which ( |4.4| ) hold is shown. 

Equation ( |4.3P is easily retrieved. Clearly, in this case 7 = a + ^ e (1, 2] thus furnishing the system 
with a superdiffusive regime. Furthermore the related coupled fractional system is ( |3.2P where ( |4.4| ) 
hold. 



Remark 4.1. If conditions ( [4.4| l hold, the equation related to the velocity field arising from p.2P can be 
determined as 

d;'dfwXz,t) = cfd,,wXz,t), t>0, zeR. (4.5) 

Here the semigroup property is no longer valid. 



Remark 4.2. Considerations similar to those in Remark 4.1 in the case of 



^e(l,2), a e (0,2-/3], (4.6) 

lead to the following results. The semigroup validity region is pictured in Fig. [I] region C. In this case we 
obtain 

d1w'{z, t) = cfd,y[z, t), t > 0, z e R, (4.7) 

and 

dfd;'pXz,t) = cfd,,pXz,t) t>0,zeR. (4.8) 

Concluding we note that the cases for which the semigroup holds, are not exactly a generalisation of 
the classical case (a = 1, ^ = 1). However we will give in the last section a physical interpretation of 
these equations in relation to the original framework, i.e. equations of velocity and mass conservation 
of one-dimensional compressible fluid. 



Remark 4.3. Here we remark that solutions to both uncoupled equations ( |4.3| ) and ( |4.5| ) (or the related 
( |4.7| ) and ( |4.8| )j can be determined by first solving the subdijfusion type equation and then by arriving at a 
solution of the sequential derivative type equation by means of the relations of the coupled system ( |3.2p . 
Indeed we have for example that the solution to the Cauchy problem 

djp'iz, t) = cfd,,pXz, t), t > 0, z e R, 

p'(z,0) = 5(z), (4.9) 



5,p'(z, 01, =0 = 0, 



is well-known and reads 



where 



P'f^'f)=^^^-r/2,i-r/2| -^ h t>0,z^R, (4.10) 



^>c,,(y) = y; ^„,l^^V K>-l,y^C (4.11) 



^ rWiKr + T])' 




2 alpha 

Figure 1: The four characteristic regions in the a^ -plane. 



is the Wright function [8. page 54]. In turn, with the initial condition w'{z, 0) = f{z), we can write 

^2 



dfw'iz,t) = -^d,pXz,t) 

Po 

^ 4dfwXz, t) = -^J^d^pXz, t) 

Po 

^ w'{z,t)=f{z)-^4d,p'{z,t\ 
Po 



(4.12) 



5 Coupling for PDEs involving sequential fractional derivatives 

In this section we present a possible way of addressing the problem of solving Cauchy problems involving 
sequential fractional derivatives. We consider the region of the a/3 -plane defined by the constraint 
(a, ^) e B (see Fig. [l]) . We seek to solve the following initial value problem: 



a°3/7(z,t) = Aa,j(z,t), t>o,z< 

/(z,0) = g(z), 

af/(z,t)|,=o = g(z). 



(5.1) 



We underline that we have requested specific initial conditions involving fractional derivatives. This 
choice will be clear in the following discussion. Here we recall that this kind of initial conditions has 
not a broadly accepted physical meaning, but some suggestions on its interpretation are already present 
in the literature (see [7]). In this section however, we are more interested in outlining a possible way to 
solve partial differential equations with fractional sequential equations more than in defining a precise 
physical meaning. 



In order to solve ( |5.1| ), we consider an auxiliary function ip{z, t), solution to 



d^dyiz, t) = A5,,(^(z, t), t > 0, z e K, 

(^(z,0) = h(z), (5.2) 

By means of ( |5.1| ) and ( |5.2| ), we can construct the coupled system 

\5«^(z,o + 'ca,/(z,o = o, 

with the initial conditions /(z,0) = g(z) and t^(z, 0) = ?i(z). Notice that functions g^z') and h(z) are 
linked together. Indeed from ( |5.3| ) we have that 



A 



3f/(z,OL=o = — 3,c^(z,0) 



^ g(z) ^ 



A 



-3,(^(z,0). 



By inserting the initial condition of ( |5.2| ) into ( |5.4| l we immediately obtain 

A 
g(z) = D,h{,z). 

K 

Similarly, from ( |5.3P , 

3^Xz,OL=o = -'c3,/(z,0). 
Considering ( |5.2| ) and ( |5.1| ) we arrive at 

h(z) = -/c3,/(z,0) = -/cD,g(z). 
Summarising, the coupled Cauchy problems becomes 

^3,Vfet)+^5,^(z,t) = 0, 
^ 5«i^(z,0 + k3J(z,0 = 0, 
VfeO) = g(z), -g{z) = -^^Bji{,^), 
c^(z,0) = h(z), h(z) = -KD,g(z), 

t > 0, z G R. Now, the application of the Fourier transform 



e(",t)^ 



e'"^CU,t)dx 



to (15. 8D leads to 



I 3"i^(6j, t) = -/ciw/Cw, t). 



/(w,0) = g(«), |(6;) = -^i6;h(«), 
(^(w, 0) = h((D), h{^oj) = —Kicog[co). 



In turn, by passing to the Laplace transform 

C(X,S): 



e-^'CU,t)dt 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



we have 



\s"(p[co,s) — s"~^/i(w) = —Kicof[cL),s). 

By substitution we can write that 

7 - 1 

s^f(,co,s)-sl^-^g(,co) = ico (s^-^Kco) - Kicof(,co,s)) — 

K 



(5.12) 



(5.13) 



fromJslffl 



^a+p-1 



s"+'7(w,s) -s"+P-^g(w) = s^-'gf^co) - Xco'-f{cc,,s) 



^ /(w,s) = g(w)- 



a+/5-l 



+ g(w)- 



We now make use of the following Laplace transform: (see Mathai and Haubold IIIOII . formula (2.2.26)) 



' ' si — n 



where 



^r,,r(y) = 1] ■ 



^rCrjr + r) 

is the Mittag-Leffler function. Therefore, by inverting the Laplace transform we obtain 
f{co, t) = g(6j)£„+^,i(-Aco2t°+^) + Kco)tPE,+p^p^,^-Xcoh'^+n 



(5.14) 



(5.15) 



(5.16) 



The inversion of the Fourier transform is now only a matter of an application of relations (14.7) and 
(14.8), page 27 of Haubold et al. O which leads to an explicit formula involving Fox functions. 



1,0 



/(z,t) = g(z)* -H,^-;' 



Xf'+P 



(l,a + /3) 
(1,2) 



(5.17) 



+ ''««' >u 



?it"+p 



(l + /3,a + ^) 
(1,2) 



t > 0, z ' 



where the symbol * represents the Fourier convolution. From the definition of the Fox function (see e.g. 
Kilbas et al. [HI) we can write 



/(z,0 = g(z) * 

+ t^g(z) * 

= g(z) * 



1 1 



r(i + 21?) 



M 



z\ 2ni J^ r(l + (a + ^)'&) \^ At«+^ 

^ 1 1 r r(i + 2'&) 

z| 27ri J^ r(l + ^ + (a + /3)i?) V At«+'' 

^ 1 If r(l + 2i?) 

/, „^fl 9.71-7 I rn 



(5.18) 






+ t'^g(z) * 



,JJ^ 2ni J^ r(l + (a + /3)i?) \^ ^Xt^. 

^1 1 r r(i + 2i?) ( 
I . .„_,_« 9.7T7 I rn 4 



d-ff 



-2i?-l 



y^^^27TiJ^r(l + ^ + (a + ^)'&) \^^^^. 



df) 



W* \ ^ ^ ( r(T) f \z\ Y\^ 



+ tPgiz) * 



ft 1 r r(T) 



V- 



2v/A^27riJ^r(l + ^ + (a + ^)^) V a/at^. 



dT 



Furthermore, by means of the Mellin-Barnes representation of the Wright function []8l formula (1.11.3)] 
we finally obtain 

f[z,t) = giz)* / W_^^_^ -^ (5.19) 



1 f Iz| 

W_a+g „_a+g ^ 

,2VAt«+^ 2 . ' 2 \^ /It^. 



+ t^g(z)* I ^ _ VV_2±g,i+^_2±g I ; : I I , t>0,Z' 



Remark 5.1. Note that ( |5.19P when g[z) = concides with the solution to the superdijfusive fractional 
diffusion equation. The reader can consult Kilb as et al. [8], Section 6.1.2. In this case ( |5.19| ) shows abo 
that the semigroup property is satisfied. 

Remark 5.2. We note that in the case (a,^) e D (or correspondingly, [a, (3) e C), i.e. when the problems 

'd^d^f{z,t) = xd,j{z,t\ t>o,z> 

/(z,0) = g(z), (5.20) 

dff[z,t)\,^o = g(.^l 

and 

'5,"+^(^(z,0 = A3,,(^(z,0, t>0,zeR, 
V(z,0) = h(z), 

are taken into consideration, a reasoning similar to that of this section, leads to the same solution ( |5.19| ) 
for f{z, t) and to the solution to the well-known fractional diffusion-wave equation for function ip[z, t), as 
it should be. 

6 Discussion 

Here we discuss the main result of this paper: finding diffusive processes from the fractional generalisa- 
tion of density and velocity equations of a simple one-dimensional compressible fluid. 

From a physical point of view we tried to understand the consequence of using a wide generalisation 
of the constitutive equations on sound wave propagation in order to take into account memory effects. In 
this way we want to show that fractional calculus proves to be a useful anal)^ic tool to describe a general 
class of phenomena beyond classical memoryless processes. Moreover, our aim is to build a complete 
theoretical framework to justify the empirical modification of the wave equation discussed in literature. 



Regarding this, for example, we recall the recent paper of Prieur and Holm II12II that relates non-integer 
power law absorption of sound waves to fractional equations; they justify the modification of the 
wave equation starting from the constitutive equations rather than considering empirical modifications 
due to observed attenuation and dispersion power laws. We also suggest an interpretation within the 
framework of fluid mechanics, recalling some relevant recent results on fractional modelling in this 
field. In their work Fellah et al. []4j] suggested the following constitutive equations of acoustic waves 
propagations in porous media: 






<t,m . - - <^6-i) 



where K^ is the bulkmodulus of air, is the porosity, p the density, p is the acoustic pressure, a(t) and 
^(t) are the dynamic tortuosity and compressibility of the air; clearly the symbol * denotes the time 
convolution 

(/*gXt)= I /(t-tOg(tOdt'. (6.2) 

Jo 

We refer to the original paper for an analysis of the physical model. What is most relevant is that they 
used fractional calculus to take into account time variability of a and p. The convolution with dynamical 
coefficients clearly implies a physical global delay in the signal propagation and the presence of memory. 
Also the experimental variability of dynamical coefficients suggests to use the time-fractional derivative 
in place of the convolution considered in ( |6.1| ). This fractional model for acoustic wave propagation 
was also experimentally validated in [5]. In this framework our model can be simply adapted. In fact, 
following [|4l, ( |6.1| l can be written as a coupled system of fractional equations on pressure and velocity 
fields in a one-dimensional model. The indices related to the order of differentiation are free parameters 
to be fit by means of experimental findings on the variability of the characteristic coefficients in the 
medium of propagation. Furthermore there is no need to assume a = ^ in our constitutive equations. 
Fractional calculus in this case emerges as a natural formalism from the formulation used in [41. This 
allows us to have a family of flexible processes through which we can deal in a general way with 
many different physical situations depending on the properties of the media of propagation. Finally, 
considering the results biefly reviewed here, we understand that fractional wave equations are relevant 
in cases of propagation in media characterised by power law frequency-dependent physical coefficients. 
From a mathematical point of view, we find two interesting results. First of all, when the semigroup 
property holds for the Caputo derivative, from the constitutive equations of the model we obtain 
diffusive or subdiffusive behaviours. Hence we can infer anomalous diffusive processes from fractional 
models of fluid such as that in ( |6.1D . On the other hand, we find a way to solve partial differential 
equations with fractional sequential derivatives based on analysing the related coupled system of 
equations. 
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